Calculation of the quantum conductance of silicon-doped carbon wire nanojunctions 
using the phase field matching theory and the tight-binding model 
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In the present work we provide the generalization of the phase field matching theory (PFMT) for 
the multi-scattering processes of the electronic excitations in a atomic wire nanojunctions. In the 
framework of the described theoretical approach, the model electronic dynamics of the considered 
systems are employed within the tight-binding approximation. Despite of using already determined 
tight-binding Hamiltonian parameters, we modify the Harrison parameters and show that covalent 
low-coordinated atomic wire systems can be described within the two-center approximation with 
a high accuracy, comparing to the first principle results. As a implementation of our model cal- 
culations, we determine the total electronic conductance of the silicon doped carbon atomic wires. 
The numerical analysis yields additionally the deep discussion of the propagating and evanescent 
fields as well as the full description of the transmission and reflection probabilities, familiar in the 
Landauer-Biittiker formalism. On the basis of obtained results we show that apart of the infinite 
periodic diatomic silicon-carbide atom wires, its finite implementations can exhibit a non-zero elec- 
tronic conductance. Additionally, we discuss various non-periodic arrangements of the carbon and 
silicon atoms in the scattering region. In a result we note that the conductance of these purely 
one-dimensional nanostuctures is a strong function of its structural properties. 
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I. INTRODUCTION 

In nature, carbon exist in a wide range of allotopic 
forms. Some of them like two-dimensional graphite [l|, 
[H spherical fullerenes Q or quasi-one-dimensional car- 
bon nanotubes @ exibit exceptional physical prop- 
erties and can be considered as a promising components 
for future nanodevices @. Recent discovery by Jin et 
al. turns the attention on another intriguing car- 
bon structure, namely the monatomic linear carbon wire 
(MLCW), which is the only stable structure among all 
forms of MLCW. In the experiment conducted by Jin 
et al. this purely one-dimensional carbon allotrope was 
produced by removing carbon atoms row by row directly 
from the graphene sheets. Apart of the earlier experi- 
ments on MLCW 3-0], method presented in Q is the 
first one which allows to obtain remarkably stable free- 
standing MLCW. 

Most often, one-dimensional atomic wires like MLCW 
are expected to have a various technological applications 
in nanoelectronics, nanomechanics, and nanomaterials 
e.g. as a interconnects between larger nanodevices 14 [. 
Due to this fact the electronic and transport properties 
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arc the key features of these small structures [l5|-[r7|. 
However, experimental research still does not provide es- 
sential knowledge about these properties of the MLCW, 
and only theoretical studies shed some light. 

Despite the fact that, MLCW were investigated for 
a long time from the theoretical point of view [Tcl|-[29j. 
they were obscured by the interest on other carbon struc- 
tures (i.e. carbon nanotubes) until the beginning of the 
last decade. It has been shown that, from the structural 
point of view MLCW can undergo the formation of ei- 
ther cumulene (all atoms connected by double bonds) or 
polyyne (alternating single and triple bonds) wires [3], 
fH, [H, [H, [ll|. However, there is no straightfor- 
ward answer which of these two structures is favorable 
one. Experimental studies does not give any satisfying 
answer on this problem, while the theoretical results de- 
pends on the applied computational method. The den- 
sity functional theory (DFT) calculations predicts double 
bond structures (32j , (33j whereas ab-initio Hartree-Fock 
(HF) methods favor alternating bond systems [l8j-[2l|. 
[30| . This situation arises from the fact that DFT meth- 
ods tend to underestimate bond alternation (second order 
Jahn- Teller effect), while HF calculations overestimate it 
[30| . Recently however, first-principle calculations indi- 
cated, that both structures are stable, whereas cumulene 
presents greater stability due to the higher Gibbs free 
energy value (34[. Furthermore, in the work by Zhang et 
al. the cumulene was predicted to be extremely strong 
pur ely one-dimensional material along the axial direction 

Moreover, on the basis of the first-principle calcula- 
tions the cumulene MLCW are expected to be almost 
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perfect conductors 34|-j45j, even better then a linear gold 
wires |32| . [46| . while polyyne are semiconducting 34 |. 
This feature proves the fact that the physical properties 
of a given material in one-dimension may be drastically 
different from those in the bulk case. In particular, three- 
dimensional diamond structure of carbon is a good insu- 
lator, while quasi-one-dimensional carbon nanotubes can 
be either semiconducting or metallic. It is also worth to 
notice that MLCW may exibit conductance oscillations 
referred to the even or odd number of atoms in the wire 

m, m. 

In our work we would like to turn readers' attention 
on the problem which to our knowledge remains still un- 
solved, namely the behavior of the electrons in the disor- 
dered MLCW. Our interest in this issue arises from the 
fact that structural or chemical impurities are common in 
physical systems and have a significant impact on a var- 
ious properties of nanomaterials [I?}, [H|. In particular, 
the foreign chemical defects injected into to the perfect 
structure allow to control the total conductance across 
the nano junction (49) . 

The consequences of such disorders of defects on the 
electronic properties of atomic size systems are usually 
viewed in the terms of Landaucr-Buttikcr theory. Such 
representation allows one to connect the classical scat- 
tering theory and the conductance of the system in the 
framework of transmission and reflection processes. We 
notice that several approaches have been developed in 
order to calculate the transmission and reflection cross 
sections, where the most popular based on the first- 
principle calculations [Hoj- [53| or semi-empirical methods 
using the Green's function formalism JHj], [55j]. In the 
present work we investigate the scattering processes on 
the basis of phase field matching theory (PFMT) (56| . 
[57| originally developed for scattering of phonons and 
magnons (58| . Our theoretical model rely upon the ma- 
trix calculus where the Bloch states of an ideal leads are 
matched to the scattering region wave functions. In this 
approach the electronic properties of the system are de- 
scribed in the framework of the tight-binding formalism, 
which have been previously successfully applied in the 
electronic transport calculations across various nano junc- 
tions [H| -[6l|. In particlar, in our work we employ 
the appropriate Slatcr-Koster [62j-typc Hamiltonian pa- 
rameters calculated on the basis of Harrison tight-binding 
theory (HTBT) [H|]. in this manner, the PFMT can 
be considered as a transparent and efficient mathemati- 
cal tool for calculations of both electronic and transport 
properties for a wide range of atomic sized systems. 

Presented work is organized in the following man- 
ner. In Sec. HH we give the detailed discussion of the 
model formalism which we need later on. The numerical 
results which incorporates propagating and evanescent 
electronic states for the disordered MLCW are presented 
in Sec. IIII1 Finally our considerations ends with the perti- 
nent conclusions and appropriate appendices which sup- 
plements the theoretical model presented in Sec. |TTJ 



II. THEORETICAL MODEL 

A. Model dynamics 

The disordered carbon atomic wire (schematic repre- 
sentation depicted in Fig. (H}) is described by the follow- 
ing tight-binding Hamiltonian block matrix 
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(1) 

which is defined for overall N x atoms in unit cell (UC) 
and Ni different orbitals per atom. In Eq. ([1]), Ejj are 
the matrices composed of both diagonal (e™' a ) and off- 
diagonal (h™'™^) Hamiltonian matrix elements for iV-th 
UC, while matrices Hij contains only off-diagonal ele- 
ments for interactions between different UC's. Each di- 
agonal parameter for a n-th atom in UC is characterized 
by the lower index I which takes the values of the an- 
gular momentum character and index a which denotes 
type of the atom (C or Si). The /i7i"nf P arame ters de- 
scribes the m-typc bond connections (m = cr, tt) between 
I and I' nearest-neighbour states, where index j3 refer to 
the types of interacting neighbors (C-C or Si-C). The 



h.v 



matrix elements are consistent with the Slater- 
Koster convention [62| and expressed in the framework 
of Harrison's tight-binding theory (HTBT) by [H 



h ll>m = VIM 



'm e dp 



(2) 



where rji^^ m are dimensionless Harrison's coefficients, m e 
denotes electron vaccum mass and the dp is the inter- 
atomic distance. Explicit forms of the Ejj and Hjj ma- 
trices are given in Appendix fff] 

In our calculations the single-particle electronic wave 
functions are expanded in the orthonormal basis of 
atomic wave functions <pi(r) as follows 

*(r, k) = ^c 4 (x n - XAr,k)0 ; (r - x n - x N ). (3) 

InN 

In Eq. (J3|), k is the real wave vector, x n denotes the 
position of the n-th atom in UC and Xn gives position 
of the UC in the atomic wire. For the ideal leads (see 
Fig. I|T])) the wavefunction coefficients ci(x n — xn,)s.) in 
the consecutive UC's are described, in the sense of Bloch- 
Floquct theorem, by the following phase relation 

ci(x n -x jv+i,k) = zci(x n -x at, k), (4) 

where parameter z is the phase factor of the form 
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FIG. 1: Schematic representation of the finite silicon doped 
carbon wire nanojunction between two semi-infinite quasi- 
onedimensional carbon leads. For the purpose of electronic 
transport properties we distinguish the so-called, irreducible 
region, and matching domains (please see section Ul Cl for more 
details). Moreover, the binding energies for a given atomic 
site and the coupling terms between neighbor atoms with cor- 
responding interatomic distances are depicted. Please note, 
that in order to keep figure transparent we are omitting the 
n and n' indices for the coupling parameters. 



z± 



(5) 



which corresponds to the waves propagating to the right 
(+) or to the left (-) direction. 

In practice, the assumed orthonormal minimal basis 
set allows us to neglect the atomic wave functions overlap 
integrals, and write the N x x Ni size eigenvalue problem 
for propagation waves in the following manner 



(EI-M d )c(x N ,k) =0, 



(0) 



where E stands for the electron eigenvalues, I is the 
identity matrix, while the dynamical matrix Md contains 
Hamiltonian matrix elements and z factors. The c(xn, k) 
is the N x x Ni size vector defined as 



c(x N ,k) 



( c s (xi - x N ,k) ^ 
c Pm (xi - xjv,k) 



\ ci(x n - x N ,k) J 



(7) 



In general, Eq. © gives N x x Ni eigenvalues to- 
gether with corresponding eigenvectors determines the 
electronic structure of the system. One can notice that 
the orthonormal character of the atomic orbitals may 
results in inadequate description of the considered elec- 
tronic properties. However, what can be seen later, the 
proper choice of the binding energies and nearest neigh- 
bor coupling terms allows us to to obtain a very good 
agreement with the density functional calculations. 



B. Evanescent states 

The description of a disordered atomic wire is addition- 
ally supplemented by the evanescent electronic states. 
This situation arises from the fact that defect site breaks 



the perfect periodicity of the atomic wire and alters the 
wave functions amplitudes, what prevents proper formu- 
lation of the problem only with the use of pure Bloch 
factors ([5]). 

Depending on the complexity of a given electronic 
state, evanescent waves arc defined by the phase factors 
for a purely imaginary wave vectors k = in as 



z± 



(8) 



or for complex wave vectors k = K\ + 1K2 as 

Z± = e T(i«i-«2)x„_ (g) 

The phase factors of the form (JSj) and ([9]) respectively de- 
fines exponential and oscillating evanescent states which 
both decay to the right (+) or left (-) direction. Please 
note that given Z-type evanescent wave corresponds to 
the energies beyond the energy band described by the 
propagating wave in this state. 

The functional behavior z(E), for both propagating 
and evanescent states at given TV, is described, on the 
basis of Eq. (|4]) and (|6]), by the generalized eigenvalue 
problem for z 



EI — Eat. w Hat,at_i 
I 



1 f XN > k l ]=0, 



(10) 



which is twice as big as the previous eigenvalue problem 
of Eq.©. In general, Eq. (flU]) gives 2N x Ni eigenvalues 
as a N x Ni pairs of z and Please note, that some 

of the solutions are trivial or non-physical. Additionally 
for systems with more than one atom per UC, matrices 
Ha^at-i and Hjy N-1 are singular. In order to obtain 
only non-trivial solutions, eigenvalue problem (TT]) is re- 
duced to the appropriate 27V; size problem, by using the 
partitioning technique (please see Appendix [B]). More- 
over, only solutions with \z\ = 1 (propagating waves) 
and j 2; I < 1 for Im[z] = (evanescent waves) are retained 
as a physical ones. 



C. Phase field matching theory 

Let us now consider the scattering problem at the de- 
fect domain (TV 6 [0, D — 1]) of the quasi one-dimensional 
wire nanojunction as presented in Fig. [TJ In particular, 
the representative set of the atomic sites includes the ir- 
reducible domain which interacts with its first boundary 
neighbors of the left and right leads, respectively. The 
general concern is to calculate the wave function coeffi- 
cients for each of the sites within the representative set. 

The Ci(x n — XN, k) coefficients of interest are described 
by the Eq. ([5]), however for the purpose of our calcula- 
tions it is convenient to rewrite Eq. ([S]) in the terms of 
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Eq. (JT)). Hence, more familiar form of the equation of 
motion for a given N is expressed as 

(EI - E N , N - Hjv,jv-i - Hj ViAr _ 1 )c(a;Ar, k) = 0, (11) 

Further, Eq. (fTTj) can be written for each point in the 
real space which is included in the representative set of 
the atomic sites. This procedure generates Nd + 2 equa- 
tions of motion of the form (|lljl and iVd + 4 corresponding 
and independent c(a;Ar,k) vectors. At this point, we no- 
tice that such set of equations cannot be solved directly, 
since the considered unknown wave function coefficients 
always outnumbered the equations of motion. Nonethe- 
less, within the phase field matching theory, this alge- 
braic problem can be reduced into the solvable size. 

In the first step the considered problem is treated as a 
scattering process within the Landaucr-Biittiker formal- 
ism. Assuming that the incoming wave propagates from 
left to right the wave function coefficients on the left and 
right side of the perturbed domain are written as 



c L (a;jv,k) = c ; (k)z 



-N 



V 



c v (\a)z N ri,v for N < -1, (12) 



c R (x N ,-k) =5^ci»Ck> Ar ii,,/ for N^D, (13) 



and reformulate its initial form (|14[) in the following way 

/ z 2 \ 

z 
1 



;(k) 



1 o 



V o 

c(x ,k) 
c(xi,k) 
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c(a;u_i,k) 



( c,(k) 
Q(k)^- 1 




V J 



(15) 



On the right hand side of the Eq. (|14j) , the rectangular 
symmetric sparse matrix has the (D + 4) x (D + 2) size. 
Remaining two vectors decouple the c(k) coefficients into 
the scattering and incident terms, respectively. Finally, 
vectors riji and tj,;/ are mapped onto the reflection and 
transmission parameters, as 



r M /(k) = y^c;/(k)r;, 



(16) 



where c;(k) and c;<(k) denotes the the eigenvectors of 
the dynamical matrix Md at given energy for electrons 
incoming in Estate and being reflected or transmitted 
in i'-state respectively. These eigenvectors correspond 
directly to the eigenstates of the perfect leads and are 
calculated of the basis of Eq. ©. However, in order to 
consider also the evanescent states it is required to per- 
form the computations by substituting the solutions of 
Eq. ([TDD into the Eq. ©. Finally, in the Eqs. {HJ and 
(|13p . the Tiji and tiy parameters are the components of 
Hilbcrt space and describe the reflection and transmis- 
sion scattering processes, respectively. 

Following Eq. ([12]) and Eq. (fig]), solutions of the 
Eq. (jT2"j) the set of the equations of motion can be com- 
pletely determined in the terms of the riji and tij< coeffi- 
cients. In particular, it is possible to group the unknown 
c(xn, k) vectors into the one Z? + 4 size super- vector c(k) 



c(k) 



f c(x_ 2 ,k) ^ 
c(x_i,k) 



c(x D ,k) 
V c(x D+1 ,k) / 



(14) 



*M'(k) 



V 



ci>(k)t u , 



(17) 



In a result, Eq. (|15[) reduces the q(x„ — xn, k) unknowns 
into the desirable number and allows to generate the fol- 
lowing inhomogencous set of equations 



/Mi,! Mi, 2 

M 2 ,i M 2 , 2 




/ Y i,i' 

c(x ,k) 
c(xi,k) 







c(iTD_2,k) 

c(a:£)_i,k) 











\ 



M.£>+2,D+1 Mj) +2 ,fl+2 / 








(18) 
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where block matrices Mij are composed of the Hamil- 
tonian matrix elements and corresponding phase factors. 
Moreover, vector which incorporates the M{ n and M| n 
elements, regroups the inhomogeneous terms of the in- 
cident wave. The explicit forms of the Mij and M?" 
components are presented in Appendix [Cl 

In practice, Eq. (fT5| can be solved using standard 
numerical procedures, yielding the c\(x n — a^k) coeffi- 
cients inside the defect region (N € [0,D — 1]), and 2L 
reflection ri^(E) and transmission tij>(E) energy depen- 
dent coefficients. 

We remind the fact that the transport processes in the 
Landauer-Buttiker theory are presented within the scat- 
tering matrix formulation. The elements of such matrix 
are respectively the reflection Ri^(E) and transmission 
Ti^i (E) probabilities, which can be determined on the 
basis of riji(E) and tijr(E) coefficients. However, it is re- 
quired to normalize these coefficients with the respect to 
their group velocities (vi) in order to obtain the unitarity 
of the scattering matrix. In what follows, the reflection 
and transmission probabilities are written as 

R u ,{E) = V -^\n, v {E)\ 2 , (19) 



1V(£) = — \t hl ,{E)\\ (20) 
vi 

where vi denotes group velocities of the waves incoming 
in state I, whereas the vi* are the group velocities of the 
waves reflected or transmitted in channel V . Discussed vi 
parameters of a given / state are calculated on the basis of 
the altered dynamical matrix (please see Appendix [Pi for 
the details). In addition, we note that for the evanescent 
states, group velocities are set to zero. This indicates the 
fact that evanescent waves do not contribute to the elec- 
tronic transport, nevertheless they are required for the 
complete description of the multi-scattering processes. 

Further, using expressions (|2Tj) and (|22|) the reflection 
probability for a wave incoming in a /-channel (Ri(E)) 
and the total reflection at the defect region (R(E)), re- 
spectively, are given by 

Rl(E)=Y,Rl,l'(E) and R(E) = ^R l (E). (21) 

V I 

Similarly, by summing over the output and then input 
channels the corresponding transmission probabilities, 
Ti(E) and T{E), can be written as 

T 1 (E)=Y,T 1 m(E) and T(E) = ^Ti(E). (22) 

V I 

The Ti(E) and T(E) probabilities are very important in 
the discussion of the scattering processes since they di- 
rectly correspond to the experimentally measurable ob- 
servable. Likewise, the total transmission T(E) allows 
to calculate the overall electronic conductance. In this 



work we assume the zero bias limit and write the total 
conductance in the following way 

G(E F ) = G T(E F ), (23) 

In Eq. (|23l) . Go is the conductance quantum and equals 
2e 2 /h. Due to the Fermi-Dirac distribution the G(Ep) is 
calculated at the Femi energy of the perfect leads band 
structure, since electrons at this level give the only im- 
portant contribution to the electronic conductance. Itself 
the Fermi energy can be determined using various meth- 
ods. In the present work, the Ef value is given on the 
basis of the density of states calculations. For a compu- 
tational details please see Appendix [Ej 



III. NUMERICAL RESULTS AND DISCUSSION 

A. The tight binding model and basic electronic 
properties 

In this section we present the electronic structure cal- 
culations which allow us to establish our choice of tight- 
binding parameters and check its possible applications 
to the considered wire nano junctions. The correctness 
of our results is validated by the comparison with the 
first-principles results after [32j], 

In the first step of electronic structure calculations 
within the tight binding approximation, it is desirable 
to carefully choose the correct minimal basis set of the 
atomic wave functions. We note, that usually carbon is 
characterized by the 2s and 2p valence orbitals, wheres 
silicon by the 3s and 3p atomic wave function types. Such 
scheme is assumed in our calculations and gives us four 
different orbitals types (namely s, p x , p y and p z ) for each 
type of the atom. Additionally, as we already stated in 
Sec. Ill A| the described set of atomic orbitals is assumed 
to be of orthonormal character. However, in contrast to 
our choice, the electronic and transport properties of the 
nanoscale systems are usually described in the framework 
of the sophisticated non-orthogonal tight binding meth- 
ods. In our opinion, such tight binding models may arti- 
ficially complicate the calculations in the case of the low- 
coordinated systems like mono- or di-atomic wire nano- 
junctions. On the basis of other experimental and theo- 
retical studies one can easily notice that the atomic wires 
have relatively simple structural and electronic charac- 
teristics. These facts proves, that for the calculations of 
electronic and transport properties of discussed atomic 
wires it is reasonable to consider the orthonormal rather 
then non-orthogonal tight binding model. 

In our case the electronic states of the considered nano- 
junctions arc described in the framework of the tight 
binding approach proposed by Harrison. This approach 
allow us to introduce previously described orthonormal 
simplifications in the theoretical formalism and still keep 
the physical character of the tight binding parameters. 
We would like to turn readers attention on the fact that, 
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in some models the tight binding parameters are consid- 
ered more like fitting parameters then physical observ- 
able. In our opinion, such approach is incorrect, since it 
do not provide full physical discussion of the electronic 
states. 

In its original form the HTBT model is considered as a 
transferable method for tetrahedral phases (63[. In par- 
ticular, the values provided by the classical form of the 
HTBT, results in only qualitative description of the elec- 
tronic states for the atomic wires composed of the low 
lying chemical elements (63j . However, one can put aside 
the universality of the HTBT approach and adjust these 
parameters to the experimental or first principle results 
by varying their values. 

In the present work the classical HTBT parameters 
has been rescaled in order to match the density func- 
tional theory results presented in 32 1, [64| . The values 
of the TB parameters utilized in our work along with 
its comparison to the classical values by Harrison are 
presented in Tab. |U We would like to indicate that in 
Tab. U both sets of parameters (the classical and re- 
parametrized ones) have the same relation to the phys- 
ical quantities. In particular, the binding energies have 
fixed negative values which are relative to the energy of 
the electron in isolated hydrogen atom. Additionally we 
note, that the values of the on-site Hamiltonian matrix 
elements for states p x , p y and p z are the same and equal 
£p' a . In contrast to the e™'" parameters the off-diagonal 
distance dependent elements are calculated on the basis 
of Eq.([2]) and can be either positive or negative. This sit- 
uation arises from the fact that one have to require the 
cohesive character of the interatomic forces. With the 
reference to this assumption, the interactions between 
orbitals of the same sign causes overlap of those orbitals 

to be positive, so the parameter h™'™^ must be negative. 
In analogy, negative overlap between orbitals of different 
signs provide positive value of the off-diagonal elements. 
The off-diagonal Hamiltonian matrix elements for inter- 
actions between s and p-type states are included only 
between s and p x orbitals, since interactions between s 
and p y or p z states vanishes in the one-dimensional case 
in order to the symmetry conditions. Additionally we as- 
sume the interactions between p x states as h™'™^, while 
the elements for interactions p y -p y a,ndp z -p z are equal in 
energies and assumed as /i™'™^. 

We note that, except of the p-state binding energy of 
silicon, all values of the TB parameters proposed in this 
work are smaller then in its classical form of the Harrison 
model. One can attribute this fact to the influence of 
the local environment of the low-coordinate systems, for 
which the binding energies are usually smaller then in 
a bulk case. However the value of the p-state binding 
energy of silicon is somehow greater when compared to 
the Harrison fit. In our opinion it is due to the lack 
of the additional high-energy orbitals on the basis set 
e.g. the peripheral Louie-type s* state [6^|. In the case 
of the nearest-neighbor coupling terms one can notice 
that the relations between the consecutive parameters for 



carbon and silicon arc identical as in its classical forms. 
Hence, the physical meaning of these is conserved. Again 
lower values of the re-parametrized parameters suggest 
environmental influence. 



TABLE I: The values of the tight binding parameters (in eV) 
and Harrison's fitting coefficients (dimensionlcss) proposed in 
this work and compared with the original values by Harrison 
[&|. Please note that the distance dependent h?',? param- 
eters are computed for the appropriate interatomic spacings 
dp (in A) tabulated below and assumed after [32] and [64|. 

Harrison TB parameters Our TB parameters 



a 






C 


Si 




C 


Si 








19.38 -14.79 




18.89 


-13.5 


e v 






11.07 -7.59 




10.94 


-8.38 


P 




c 


Si 


SiC 


C 


Si 


SiC 


Vs,s 


a 




-1.32 




-0.93 


-1.48 


-1.11 


Vs,p 


a 




1.42 




0.94 


1.19 


0.95 


Vp,p 


a 




2.22 




1.03 


1.18 


0.99 


Vp,p 


7T 




-0.63 




-0.59 


-0.41 


-0.62 


p 




c 


Si 


SiC 


C 


Si 


SiC 


hs,s 




-5.95 


-2.08 


-3.70 


-4.19 


-2.33 


-3.11 




a 


6.40 


2.24 


3.98 


4.23 


1.87 


2.66 


hp^p 


<T 


10.01 


3.50 


6.22 


4.64 


1.86 


2.77 


hp^p 


TV 


-2.84 


-0.99 


-1.77 


-2.66 


-0.65 


-1.74 


p 




C 




Si 






SiC 


dp 




1.3 


2.2 




1.649 



The electronic structures of the infinite MLCW, ML- 
SiW, and DLSiCW, obtained using the re-parametrized 
Harrison tight-binding terms are compared with results 
computed on the basis of classical HTBT parameters and 
previously reported DFT calculations in FigfSJ It is also 
worth to note that in the case of the first two electronic 
structures (respectively for MLCW and MLSiW wires) 
we assume only one atom per PUC, whereas the SiC 
diatomic wires consist of two atoms per PUC. More- 
over, the disused band structures are supplemented by 
the corresponding density of states results. In addition 
we note that, our tight-binding results are represented 
by the solid curves whereas the classical HTBT results 
and DFT bands (after (32[, [HI) by open triangles, and 
dots, respectively. Please note that the re-parametrized 
tight-binding electronic states arc additionally marked 
by colors due to the proper distinction between their 
different characters (such convention is assumed for all 
results in this paper). For first two structures (respec- 
tively MLCW and MLSiW nanojuntions) red and blue 
colors correspond respectively to the a and a* bands. 
These bands arises from the sp x orbital hybrids, where 
the lowest lying bands are always occupied by two elec- 
trons. Bands marked by red color has 7r character and 
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DLSiCW 



Harrison fit [55]: 
ssa 
spa 
ppa 

o ppn 
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<UA] 



FIG. 2: The nearest-neighbor tight binding coupling param- 
eters as a function of the interatomic distance. 




25 50 
DOS [a.u.] 



are doubly degenerate. Their origin lays in p y and p z or- 
bitals, allows them to hold up to four electrons. For both 
MLCW and MLSiW structures the 7r-type bands crosses 
the Fermi energy level, which means that these atomic 
wires has purely metallic character. Instead of first two 
linear wires the electronic structure of the DLSiCW con- 
sists of overall eight bands. Starting from the DALSiCC 
band structure minimum consecutive bands has their ori- 
gin in the following orbitals: C 3s (red band), Si 3s (green 
band), C 3p (blue and black bands), and Si 3p (pink and 
orange bands). Additionally for such energy levels order 
one can notice two doubly degenerate 7r-typc bands, re- 
spectively indicated by the blue and orange colors, and 
an insulating character of DLSiCW. 

We would like to turn readers' attention on the fact 
that stated metallic or insulating characters of the con- 
sidered atomic wires arc appropriate only when the wires 



FIG. 3: Electronic structures of carbon (A), silicon (B) and 
and silicon carbide (C) linear atomic wire nanojunctions, pre- 
sented over the first Brillouin zone, ka £ [— it, it]. Our results 
are represented by the solid color curves and compared with 
the first-principle results (open dots, ka £ [0, 7r]) HH, [(54l |. 
and the tight binding results obtained using the classical pa- 
rameters by Harrison (open triangles, ka £ [— it, 0]) [63| . 



are infinite. It is a well known fact that this character 
can change for the case of finite number of atoms in a 
given wire or due to the type and quality of connected 
electron reservoirs. 

Furthermore, on the basis of the Fig. [3j one can no- 
tice that, for MLCW and MLSiW, the re-paramctrized 
HTBT terms reproduces the DFT result with the quan- 
titative agreement put to the energies slightly above the 
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Fermi level. Electronic branches in the regions of high- 
energies are in qualitative agreement. In the case of DL- 
SiCW some of the electronic states, even in high-energy 
regions perfectly matches the DFT results. Our parame- 
ters in this manner constitute the most optimal set of the 
parameters for the electronic transport calculations, since 
it allows us to consider more reasonable energy ranges 
that the classical HTBT terms, and what is even more 
important, correctly reproduce the Fermi level. 



B. Evanescent states 

The electronic structure characteristics for the infinite 
MLCW, which constitute the electronic leads in our cal- 
culations, are additionally supplemented by the numer- 
ical results for the generalized eigenvalue problem (fTU|) 
in two possible representations. The top row of the 
Fig. (|3|) presents the three-dimensional representation 
of the z(E) generalized functional behavior for all elec- 
tronic states symmetries. Please note, that electronic 
states with \z\ = 1 correspond to the propagating waves, 
whereas the solutions for \z\ < 1 and Im[z] = describe 
evanescent waves. The moduli of the complex z factors 
are presented in the bottom row of the Fig. ([3]). In 
this case each plate denotes propsgating states, while the 
evanescent waves are represented by the exponential-like 
curves. We would like to emphasize the fact that results 
presented in Fig. ([3]) provide more complete view on 
the electronic properties of a given system then a typical 
band structure representation as in Fig @ , since it repre- 
sent both propagating and localized states. Such general 
representation is necessary for the proper discussion of 
the scattering problem in a disordered systems. In partic- 
ular, on the basis of Fig[J] (B) one can notice that, as we 
stated in section I, the evanescent states of given states 
are always outside the region defined by the propagat- 
ing waves in this state. However, evanescent solutions to 
the eigenvalue problem (|10p overlap the energy ranges of 
propagating states of another ^-typc character. This fact 
clearly indicate the importance of the evanescent states 
in the description of electronic transport properties for a 
nanoscale systems. 

In what follows, the energy range considered in our cal- 
culations is restricted by the electronic structure bound- 
aries. It creates a energy window of physically acceptable 
solutions. This energy range is marked by two grey ver- 
tical lines in FigfJ] (B). 



C. Transport properties 

In the last step the electronic transport properties of 
the silicon-doped carbon wire nanojunctions between two 
semi-infinite one-dimensional carbon leads have been cal- 
culated using the PFMT method. 

In Fig. ([S]), the group velocities of electrons in the 
MLCW perfect leads are presented. It is these group ve- 




-20 20 
E-£ p [eV] 

FIG. 4: Three-dimensional representation of the functional 
behavior z(E) on complex z plane (A) along with the evo- 
lution of the absolute value of the complex z factors as a 
function of the energy E for MLCW (B). Given electronic 
state symmetries are marked by different colors correspond- 
ing to the electronic structure characteristics from Fig. ([JJ. 
Please note, that the results are presented for the intervals 
E — [—50, 50] eV and Im[z] = [0, 1]) in order to receive trans- 
parent representation. 



locities of the electrons incident from the perfect leads, 
which arc used to normalize the scattering matrix. 
In particular, five different arrangements of the silicon 
atoms in the nanojunction region have been considered 
(please see Fig. ©). The first three represents peri- 
odic DLSiCW composed of respectively 1, 2, and 3 Si-C 
atomic pairs. Next case correspond to the nanojunction 
wire system composed of 3 carbon and 3 silicon atoms, 
but without the long range periodic arrangement. Fi- 
nally, fifth system incorporates 5 silicon atoms and only 
one carbon atom in the middle. 
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FIG. 5: Group velocities over the intervals of propagating 
states for MLCW perfect leads. 



Case 1 




FIG. 6: Schematic representation of the five different arrange- 
ment of the silicon and carbon atoms on the nanojunction wire 
region, which are considered in this section. All presented 
nanojuctions are assumed to be held between two semi-infinite 
quasi-onedimensional carbon wires. Here the irreducible re- 
gion for case 1 and 2 is marked by the dotted lines, whereas 
for remaining cases only irreducible regions, without lead's 
atoms, are presented. 

The calculated transmission and reflection component 
scattering cross sections for each of four available trans- 
port channels are presented in Fig. ([7]). Each row of the 
subfigures correspond to the given considered arrange- 
ment of the silicon and carbon atoms on the wire nano- 
juction i.e. case 1 - (A)-(C), case 2 (D)-(F), case 3 (G)- 
(I), case 4 (J)-(L), and case 5 (M)-(O). The red and green 
solid lines represent the transmission and reflection com- 
ponent cross sections, respectively. The unitarity condi- 
tion (light blue lines) is used systematically as a check on 
the numerical results. Finally, the Fermi level is marked 
by dashed lines and set as reference at zero energy. Once 
again we impose the fact that, on the zero bias limit, the 
total conductance is read at this Fermi level. 

On the basis of obtained results we observe that, the 



transmission spectra presents strongly localized scatter- 
ing resonances, which number differ depending on the 
number of silicon atoms on the nanowire. In what fol- 
lows, valence cr-state (state 1) exhibit small transmission 
in each of the five considered nanojunctions. Only in the 
cases 2 and 3 narrow transmission peaks can be noticed. 
At the zero bias limit the most significant channels are 
constituted by the degenerate 7r-type states (state 2 and 
3), since only these states crosses the Fermi level. Apart 
of this fact, also the characteristics of the high-energy 
states (state 4) may prove to be useful in the discussion 
of the electronic transport across the silicon-doped car- 
bon wires. In each of the five cases, the transmission 
around the energy of 5 eV is close to the unit. This fact 
will be important, when applied bias rise the energy of 
the incident electrons, allows them to be transmitted at 
higher, then the Fermi level, energies. 

The total conductance of the considered systems has 
been depicted in Fig. flS), as a function of energy (red 
solid lines). Wc note that G(E) is expressed in the Go 
units of the conductance quantum. Furthermore, we have 
plotted the conductance of the infinite MLCW (solid light 
blue line) and the Felmi level position (dashed lines). 
The total conductance spectra directly correspond to the 
transmission probabilities of the appropriate nanojunc- 
tion wires presented in Fig. ([7]). However, since 7r-type 
states (states 2 and 3) are degenerate, the value of the 
G(E) in the energy range of these states is in fact, the 
summation over these two identical channels. In addi- 
tion, we notice that, the most significant contribution to 
the total conductance in the zero bias limit is given by 
these degenerate 7r-type states, since only these states 
of the perfect leads crosses the Fermi level. In partic- 
ular, for the three first cases which represent the DL- 
SiCW of different lengths, the total conductance at the 
Fermi level is non-zero. This observation is in contradic- 
tion to the insulating character of the infinite DLSiCW. 
On can connect the non-zero value of the conductance 
of a short DLSiCW nanojunction wires to the band gap 
around the Fermi level in the infinite DLSiCW systems 
(please see Fig. [3] (C)). Further, this band gap is re- 
lated to the difference between the binding energies of 
the silicon and carbon atoms, respectively. Mentioned 
energy difference creates the effective potential barrier 
for the electron crossing the nanojuction region from left 
to the right lead. When the wire length is increased the 
incident electrons encounter periodic repetition of such 
silicon-carbon potential barriers, and the cumulative ef- 
fect of the conductance value decrease when the length of 
the DLSiCW nanojuntion increases, is observed. In our 
results this effect can be observed both in figures ([7]) and 
([8]). where the transmission and consecutively total con- 
ductance spectra around the Fermi level exhibits minima. 
These minima are stronger when the wire nanojunctions 
are longer. 

The two remaining silicon-carbon nanojuctions of cases 
4 and 5, respectively, are no longer of the periodic types. 
In the case 4, the number of the silicon atoms on the 
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FIG. 7: The transmission and reflection probabilities arcross five considered types of the silicon-doped carbon wires between 
two semi- infinite one-dimensional carbon leads. The arrangement of the subfigures is following: (A)-(C) - case f , (D)-(F) - case 
2, (G)-(I) - case 3, (J)-(L) - case 4, and (M)-(O) - case 5. The Fermi level is set at the zero energy reference position. 



nanowire is the same as in the case 3, however there is no 
more long range order in the arrangement of the carbon 
and silicon atoms. Such structural defect, for relatively 
short wires, completely decreases value of the G(E) at 
the Femi level to zero. Further, in the fifth case, the 
number of the silicon atoms outnumber carbon atoms in 



the nanojunction region. Similarly as in the case of the 4, 
the total conductance at the Fermi level is zero, however 
the overall shape to the G(E) spectra is different. I this 
case the zero value of the conductance is imposed by the 
large number of the chemical defects constituted by the 
silicon atoms on the nanojunction wire. 
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FIG. 8: Total electronic conductance spectra as a function 
of energy for silicon doped carbon wires (red solid lines). 
Furthermore, the conductance spectra of the perfect infinite 
MLCW nanojunction is depicted (light solid lines). Each 
G(E) function is expressed in the Go units, and the Fermi 
level is set at zero energy reference position (dashed line). 



IV. CONCLUSIONS 

In the present work, the electronic and transport prop- 
erties of the silicon-doped carbon wire nanojunctions 
wires arc discussed within the phase field matching the- 
ory on the level of the tight-binding approximation. The 
minimal basis of the atomic wave functions is assumed 
to be composed of four different atomic orbitals, namely 
the s, p x , p y , and p z types. 

In the first step of our discussion, we present new val- 
ues of the tight-binding parameters for the prefect infinite 
carbon (MLCW), silicon (MLSiW), and silicon-carbide 
(DLSiCW) nanojunction wires, by matching these to the 
previously reported density functional results. During 
the fitting procedure we benefit form the Harrion tight- 
bidning scheme, which allows us to neglect the overlap 
Hamiltonian terms, and assume the two-center approxi- 
mation for the Hamiltonian interaction integrals. 

Such set of the tight-binding parameters is then used 
in order to calculate the electronic transport properties 
of the silicon-doped carbon wire nanojunctions wires. In 
particular, we consider five different cases of the atoms 
arrangement on the nanojunction region. We show, that 
apart of the infinite insulating DLSiCW nanojunctions 
its finite implementation exhibit a non-zero conductance 
values. This fact is explained by arguing that the energy 
difference between the binding energies of the silicon and 
carbon atoms results in the effective potential barrier for 
the electron crossing the nanojuction region from left to 
the right lead. Next we consider the structural defect in 
the diatomic silicon-carbon arrangement of the atoms on 
the nanojunction wire. By changing the position of two 
silicon and carbon atoms we show that the total conduc- 
tance of the systems drops when comparing to the corre- 
sponding periodic configuration of the atoms, resulting in 
the zero value of the total conductance. In the last step, 
we present the case when the silicon atoms outnumber 
the carbon atoms on the nanojunction region. Again, 
the zero conductance of the system is observed. This is 
due to the chemical defect influence of the relatively large 
number of the silicon atoms. 

In summary, the PFMT method prove to be useful tool 
for the quantum transport calculations, and may be con- 
sidered as a alternative to the non-equlibrium Green's 
function methods. Furthermore, it is especially inter- 
esting for the treatment of complex systems presenting 
multi-channel conductance, and for the calculation of the 
conductance properties of nanojuctions which can exhibit 
electron-phonon interactions with or without chemical 
and structural defects. 
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Appendix A: Explicit forms of the E, j and Hjj 
matrices 



On-diagonal (Ei,j) and off-diagonal (H, i3 -) Hamilto- 
nian matrices are in agreement with the following def- 
inition of the one-dimensional tight-binding Hamiltonian 



Both Eq. (|A2[) and (|A3[) denotes N x Ni square matrices, 
where matrix (|A3[) is upper triangular. In this manner 
component matrices (|A4[) and (|A5[) are of the dimension 
iV; x Ni. Additionally, matrix £j'j' always denotes diag- 
onal matrix, while hi>ji matrix is much more complex, 
with possible non-zero elements at every position. Please 

note, that some of the h 7 ?'™^ elements can vanish due to 
the symmetry conditions and simplifies the notation of 
the h,v i> matrix. 



H 



22 e ?' a ^l,l'^n,n^N,0 



InN 



+ tiij?^ (Sn,n±l$N,0 + ^n,n±l^JV,JV±l)) (Al) 

where symbol <5 denotes Kronccker delta. Analogically to 
the Eq. (TTJ) , Hamiltonian (|A1[) is defined for general case 
of overall N x atoms in unit cell and Ni different orbitals 
per atom. After Eq. (|A1[) . the general explicit form of 
Ejj matrix is given by 



h 2 ,i £2 



E; 



V h„,i • ■ ■ 
while matrix Hij by 

/0 ••• hi >; 



H 









\o 



where 



H',j' 



\ 



hn,n— 1 



hi >n _i hi „ ^ 

h2,n-l h2,„ 



/ ^ 



£ n ' a 

Pa; 



/ 



\ 



-J-l u 

V / 



(A2) 



(A3) 



(A4) 



Appendix B: Partitioning technique 

The partitioning technique is suitable method which 
allows to avoid singularity problem of the Hjv,jv-i and 
H]y j V _ 1 matrices and calculate only non-trivial solutions 



of eigenvalue problem (|10p . Detail discussion of the par- 
titioning technique is presented in [o5| and this section 
gives only our short remarks on this method. 

Following studies from [o5j problem (|10[) is partitioned 
into two parts, of respectively D\ — Z? 2 and -D 2 sizes, 
where 



and 



Di = N X N U 



D 2 =N n N h 



(Bl) 



(B2) 



In Eq. (|B2[) . parameter N n stands for the order of 
nearest-neighbour interactions assumed in calculations 
e.g. N n = 1 for the first nearest-neighbour interactions, 
on the basis of Eq. (|B1I) and (|B2|) . reduced 2Ni eigen- 
value problem is written as 



Ai.i 
I2.2 



Ai, 2 




B1.1 




B1.2 
I2.2 



c 2 (x N ,k) \ =Q 
c 2 {x N -i,k) 



(B3) 



At this point we correct the misprint from [€351 ] and write 
the submatrices of Eq. (|B3|) in the following form 



and 



Ai,i — SIa,2 — E 2 ,2 — E 2 ,i (-ETi.i — Ei,i) 1 Ei j2 , (B4) 



Ai t 2 = — H2,2 — E 2j i (-ETi,i — E^i) 1 Hi j2 , (B5) 
Bi,i = H22 + Hi 2 (EI 1A - Ex,!)" 1 Ei, 2 , (B6) 



/ Ln,n',/3 j-.7i.ri' ,B 

,L S,S,<7 U S,p x ,<J 

Ln,n',/3 Ln,n',/3 



jn,n' ,/3 



fy-lj'-l,m n i-lM,rn 



rn.n' ,8 
h l,V-l,r 



in,n' ,f3 
n U'm 



(A5) 



Bi 2 — Hi j2 (£li,i — Ei,i) Hi, 



(B7) 



Please note, that reduced problem (|B3[) gives 2A; eigen- 
values with 2iVz corresponding eigenvectors, what is N x 
times less then can be expect from the physical point of 
view. Nevertheless, those solutions can be easily sepa- 
rated into N x Ni eigenvalues and N x Ni eigenvectors of a 
purely physical character. 
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Appendix C: Explicit forms of the M.ij, Ml™, and 
M™ components 



The sub- matrices of the equation (fl8|). for a given i 
and j indices, are given as 



? = J 

Mi j = Ei for <( N > i > 1 
N > j > 1 




Mi 




(Cl) 



(C2) 



(C3) 



except of the sub-matrices which describe the boundary 
atoms of the system and are expressed in the following 
manner 

Mi,i = H t _ N+3 _ N+2 c P (k)z, 
M 2>1 = -tt_ N+4 ^ N+3Cl ,(k)z N - 3 , 

Mjy— i,iv = -U\ f _ itN _ 6 cy(k)z N - 3 , 

M N . N = Hjv-3,JV-4C/'(k)z. (C4) 



Finally, the M™ and M™ vector components are written 
as 



H 



t 

-N+3.-N+2 



Ci(k)z 



N-4 



and 



Mo 



-H_ W+ 4 i _JV+3Ci(k)z' 



-N+3 



(C5) 
(C6) 



Appendix D: Group velocities 



As specified in section III Cl the group velocities for 
individual states, can be calculated on the basis of eigen- 
value problem (j6j rewritten in the following manner 



[vI-V]v{x N ,k) = 0, 



(Dl) 



where v denotes the eigenvalues of Eq. (|D1[) which yields 
all required electron group velocities for each propagating 
state. Further V is the N x x Ni size matrix of the form 



V 



<9H 



(D2) 



Finally v(xat, k) stands for eigenvectors of problem (|D1|) . 
We note that usually equation (|D2j) includes the constant 
part dp/h, where h is the Planck constant. However, for 
the purpose of electronic transport calculations within 
the PFMT approach, this term can be omitted due to 
the fact that only the ratios of the given group velocities 
are important (please see Eq. ([21]) and Eq. ((22])). 
Appendix E: Density of states and Fermi energy 



In order to calculate the density of states for a given 
system we define the retarded matrix Green's function of 
the form 



G(E,k) = p + Ol-H]" 



(El) 



using the typical mathematical trick i.e. adding small 
imaginary part i£ to the energy E. Next, the total den- 
sity of states can be calculated on the basis of the follow- 
ing relation 



DOS(£0 = / Im[Tr[G(S,k)]] 



(E2) 



The partial DOS is then simply calculated on the basis of 
general equations (|Elj) . and (|E3j) . for a given eigenvalue. 



In what follows, the Fermi energy for a given system 
is determined on the basis of carrier concentration calcu- 
lations. Assuming total N e electrons per atom, the total 
carrier concentration 



C = 



Im[Tr[G(£,k)]], 



(E3) 



is computed in the energy range from the band minimum 
Emin to the band maximum E max . Similarly to density 
one states Eq. (|E3[) serve also as a basis for the calcula- 
tions of partial carrier concentration values for a given 
electronic state. 

Now, the carrier concentration of the lowest state oc- 
cupied by the N e > atoms can be calculated, to be next 
multiplied by the number N e /N e >, and give the Ct car- 
rier concentration of the total N e electrons per atom. 
Finally if C in Eq. (|E3[) is set to be equal to Ct the 
Emax is obviously the Fermi energy. 
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